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1 Introduction 

The study of charmonium decays has seen in the 
last years a remarkable progress and a renewed 
interest. New experimental measurements, main- 
ly coming from Belle, BES, CLEO and E835 
PQ have improved existing data on several ex- 
clusive hadronic decay channels leading to bet- 
ter determinations of the charmonium branch- 
ing fractions. On the theory side several heavy 
quarkonium decay observables can be studied 
using effective field theories of QCD as Non- 
Relativistic QCD (NRQCD) and its extensions 
[2]. The latter is also a recent review in general 
of quarkonium physics. 

In this letter we focus on the OZI suppressed 
decays of \c0 ■> Xc2 into two light pseudoscalar 
mesons, irir, KK, rjr) (xcO,Xc2 — ► PP) and we 
present the arguments why we have no similar 
results for J/tp, ip{nS), Xci decays to the same 
final state. We explore the possibility to calcu- 
late the contributions to the decay amplitudes 
due to the so-called chiral logarithms. These 
have the form m 2 log (m 2 ), where m is the mass 
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of a light mesorQ and these are potentially the 
largest contribution from the light quark masses. 
Surprisingly, we find that for XcO,Xc2 — > PP 
these contributions vanish. This is not necessar- 
ily the case for scalar quantities at high q 2 as 
we show with the case of the scalar formf actor. 

Terms of the type m 2 log(m 2 ) and other non- 
analytic dependence on the input parameters 
can be produced by soft meson loops [3]. The 
modern version of this method is Chiral Per- 
turbation Theory (ChPT) gUSUE], see also the 
lectures 0[S]. All of these applications require 
the octet of pseudoscalar mesons to have soft 
momenta which is not the case in charmonium 
decays. However, it has been argued that even 
for cases with pseudoscalar mesons at large mo- 
menta there are still predictions possible. This 
was first argued for Ku, decays in [5] and later 
argued to be more general and applied to K — > 
7T7T [TU] and B — ► D,n,K,r] vector formfactors 
and the pion and kaon electromagnetic form- 
factors [111112) . The underlying arguments were 
tested at two-loop level for the pion vector and 
scalar formfactor in [T2] . We refer to this method 
as hard pion ChPT (HPChPT). 



1 In the remainder of this letter m stands for m n , 
rriK and m^. 
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The underlying argument in HPChPT is that 
the chiral logarithms are coming from the soft 
part of all diagrams when we envisage a the- 
ory with hadrons with an effective Lagrangian 
to all orders. This should be able to describe all 
effects according to Weinberg's folklore theorem 
[3]. The hard part of these diagrams around a 
particular kinematical situation are describable 
by a tree level Lagrangian that is analytic in the 
soft physics. Loops using the latter Lagrangian 
thus reproduce the nonanalytic dependence of 
the underlying loop diagrams in the full theory. 
The second part is then proving that we use a 
sufficiently complete tree level Lagrangian. This 
has been done in all the previous works [51 11 01 
ITT1IT2"] by showing that higher order operators 
have the same matrix element as the lowest or- 
der operator up to terms analytic in the light 
pseudoscalar masses. This allowed the determi- 
nation of the chiral logarithms for the processes 
mentioned earlier. The same arguments are ap- 
plicable to the present decays. 

In Section [2] we define the notation and give 
the lowest order Lagrangians involving charmo- 
nium fields. Section|3]describes the relevant loop 
calculations. We note that our zero result for the 
chiral logarithm is in fact valid for all fully chiral 
singlet states of spin and 2. One consequence 
of our result is that light quark mass corrections 
in the \c0, Xc2 PP decays are not enhanced. 
We compare this statement with the available 
experimental results in Sectional 



2 Formalism 



First we shortly summarize the formalism of 
ChPT in its three-flavour version [6 . Introduc- 
tions to ChPT can be found in [7,8 . Hereafter 
we will use the same notation as in |13j . In 
ChPT it is assumed that the spontaneous sym- 
metry breaking of chiral symmetry takes place. 
In group theory language it has the pattern 
SU(3) R x SU{3) l /SU(3)r -> SU(3) V . The os- 
cillations around the vacuum are described by 
the field u £ SU{3) 



doscalar mesons 



u = exp 



y/2F 



(1) 



where is an hermitian matrix containing the 
pseudo Goldstone bosons, i.e. the light pseu- 



V 



7T 

K~ 



K° 



75' 



K+ 
K° 



(2) 



The lowest order Lagrangian describing the strong 
interactions of the mesons must satisfy the same 
symmetries of QCD and reads 



with 



C$ = ^({u^) + (X + )), (3) 



X = 2B (s + ip). 



*n and r M = 



The fields s, p, = 
are the standard external scalar, pseudoscalar, 
left- and right-handed vector fields introduced 
by Gasser and Leutwyler [5US]. s includes the 
light quark mass matrix. 

The field u and it M transform under a chiral 
transformation g^ x g R £ SU(3)l x SU(3)r as 



9ru 



}v = hug 



Li 



hu^h/ . 



(4) 



In Q h depends on u, g^ and gu and is the 
so called compensator field. The notation (X) 
stands for trace over up, down, strange. 

As anticipated in Section [T] the use of ChPT 
is limited to those processes where the pseu- 
doscalar mesons are not very energetic. Indeed 
the power counting is based on an expansion in 
m 2 / A 2 and in p 2 / A 2 , where A is around 1 GeV 
and is the scale up to which ChPT is believed 
to work. However sometimes the dependence on 
the light pseudoscalar masses outside such en- 
ergy regime is needed. In such cases the exten- 
sion of ChPT to hard pion ChPT is necessary. 
For an extended explanation of the arguments 
leading to hard pion ChPT we refer the inter- 
ested reader to [51 HU1HT] . Here we only stress the 
key features of such extension. When calculating 
the amplitudes of processes involving hard pseu- 
doscalar mesons the Feynman diagrams contain 
both hard and soft lines. The application of hard 
pion ChPT involves a separation of the first 
ones from the latter. The hard lines get con- 
tracted into a vertex of another effective La- 
grangian. The couplings appearing there must 
therefore depend on the hard quantities. The 
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soft lines instead encode the dependence on the 
soft parts, like the light masses of the pseu- 
doscalar mesons, and are used to calculate the 
chiral logarithms. Such calculation is done us- 
ing a Lagrangian built up in the same fashion 
as the one of standard ChPT since the hard 
parts must satisfy the constraints from chiral 
symmetry. However one must keep in mind that 
now the expansion that survives is in the small 
parameter m 2 /A 2 , while the one in momenta 
must be dropped. This means that adding ex- 
tra derivatives to the Lagrangian is in principle 
allowed since they are not suppressed by extra 
powers of momentum. Fortunately it often turns 
out that the matrix elements of operators con- 
taining these extra derivatives are proportional 
to those of the lowest order operator up to terms 
either analytic in m 2 or higher order. The co- 
efficient of m 2 log m 2 is thus predictable. It is 
important then that this happens in the case 
under study as well. 

The arguments in !). 10 are for a heavy to 
two light decay and the exact same arguments 
apply to the present case of charmonium to two 
light pseudoscalars. 

Charmonium states are full chiral singlets, 
i.e. they do not transform at all under the chi- 
ral SU(3)l x SU(3)r. Their kinetic terms are 
thus the standard kinetic terms for fields of the 
given spin without couplings to pions. The re- 
sults obtained here are thus valid for all full chi- 
ral singlets. 

We only study decays to two light pseudo- 
scalar mesons, in principle we could study de- 
cays to three light mesons as well but the ex- 
tension is not totally straightforward since the 
HPChPT should be applied to each kinematic 
configuration separately. Also for more particles 
the operator structure is typically much more 
involved. 

We could as well study J/tp,tp(nS) and Xd 
decays into two light pseudoscalar mesons. Again 
the reason why we have not done this is that 
it would be impossible to compare J/tp — > tttt 
with respect to J/i/j — > KK and similar for the 
others. These two decays are caused by differ- 
ent operators. The first one violates G-parity0 
it thus proceeds either through the electromag- 
netic interaction or the quark mass difference 
nid — m u . The decays to a pair of kaons have 
the same problem but instead with U or y-spin 

2 Or more general, an L = 1 state of two pions is 
in an / = 1 state. 



rather than isospin. These decays thus proceed 
electromagnetically or through the mass differ- 
ences m s — m u or m s — m^. Notice that this is 
also the reason why the decay of J/ip into tttt is 
suppressed compared to the one into KK. It fol- 
lows that a simple comparison, as the one done 
in Section [3] for x c o,Xc2, is not possible. 

We describe the XcO state with a chiral sin- 
glet, scalar field xo an d the Xc2 with a chiral 
singlet symmetric traceless tensor field X2^i/ sat- 
isfying X2/J.V = X2^, rfxitw = and on-shell 

VxX2p,v = 0. 

We remind the reader that we are interested 
in calculating the dependence of the amplitudes 
on m 2 . To predict this we will stop at the chi- 
ral logarithm level so we only calculate contri- 
butions like to 2 log m 2 . Terms of order m 2 are 
instead of higher order and thus neglected. This 
means we can neglect all effects of x± except 
for the contribution to the light meson masses 
and their effects in loop diagrams. In particular, 
effects of the light quark masses on the charmo- 
nium mass are linear in the light quark masses 
or higher order. 

It turns out that there is only one lowest 
order operator for each case 

C Xc = EiFq xo (u»u^)+E 2 F 2 X 2 • (5) 

We have added a factor of Fq, the chiral limit 
pion decay constant, to have the lowest order 
independent of Fo. We have actually added a 
few more terms for the scalar case as an explicit 
check on the HPChPT arguments. These are 

C E Xc = E 3 F 2 xl (u^)+E 4 F 2 xo ( WV M ^) . 

(6) 



3 The Calculation 

The tree-level diagrams and the loop-diagrams 
that contribute to Xco > Xc2 —> PP are shown in 
Figure [TJ To these diagrams we also need to add 
the wavefunction renormalization for the pseu- 
doscalar meson. The wave function renormaliza- 
tion for the charmonium states has no contri- 
butions of order m 2 log m 2 . The tree level result 
from the diagram in Figure [TJl) reads 

iA( Xc0 -> PP) = -{m 2 x - 2m 2 P ) x 

(2E 1 +E 4 (m 2 x -2m 2 P )) , 
iA( X c2 -> PP) = -4E 2 p^ 2 T Xflv . (7) 
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Fig. 1. The tree-level (1) and the loop-diagrams 
(2,3) contributing to the decays XcO,Xc2 — > PP up 
to the order m 2 logm 2 . The wiggly lines correspond 
to a \c while the continuous line are tt, K or r\. A 
round vertex correspond to an interaction from £ Xc 
while a box is a vertex from L-n. Notice that the 
diagrams of type (3) do not contribute with chiral 
logarithms to the amplitudes. 




Fig. 2. An example diagram involving a propa- 
gator of the heavy particle in the loop. These do 
not contribute to the chiral logarithms of order 
m 2 log m 2 . 

The amplitude is the same for the final states 
" ,K+K~ ,K°K and 7777 when the meson 
mass mp is chosen to be the appropriate one. 
T XfJ ,u is the polarization vector of \ci and p\ and 
P2 are the fourmomcnta of the two mesons in the 
final state. Here we also see that the example of 
a higher derivative term containing indeed 
gives a result proportional to the lowest term 
up to corrections of order m 2 . 

When we go over the powercounting argu- 
ments, only soft pseudoscalar propagators where 
no derivatives act on the pseudoscalar fields as- 
sociated with that propagator can give contri- 
butions of order to 2 logm 2 . As a consequence 
only pseudoscalar meson wave function renor- 
malization and the diagram in Figure QJ2) can 
contribute terms of order m 2 log m 2 and we find 
that these two contributions exactly cancel for 
all decays considered. We indeed find that the 
diagram in Figure[lj3) does not contribute terms 
of order m 2 log m 2 . 

In baryon ChPT, and other extensions with 
heavy particles, one can also get logarithmic 
contributions from diagrams involving propaga- 
tors of the heavy particle. The powercounting 
argument together with the fact that all inter- 
actions terms have derivatives show that these 
should not be present here. We have explicitly 
confirmed this with the diagram of Figure [H 
Putting in the vertices with E\ and E3 we in- 
deed find no contributions of order m 2 log m 2 . 



We conclude this part with a final remark. 
The above discussion is valid assuming that other 
charmonium states are far enough in mass, so 
that such states cannot appear as virtual par- 
ticles in the loops, i.e. their contributions are 
sufficiently hard that they can be described by 
the same tree level Lagrangian ([5]). 

We were somewhat surprised to find that 
there were no chiral logarithms of leading or- 
der in these decays, and more, since the only 
assumption that goes in is the HPChPT and 
the fact that we have a chiral singlet field. The 
same arguments go through for the energy mo- 
mentum tensor T^ v which is a spin two chiral 
singlet. The expressions for the matrix element 
of the energy momentum tensor between pseu- 
doscalar states (P\T^ V \P) are known to one- 
loop order in ChPT 14 . Our result should pre- 
dict the m 2 log m 2 parts of their result in the 
limit q 2 > m 2 . Expanding (29) and (33) in [H] 
in m 2 / 'q 2 we indeed find that there are no loga- 
rithms of order m 2 log m 2 appearing there. 

In earlier work we have found many instances 
of nonzero chiral logarithms in HPChPT. Most 
of these were for the vector formfactor but we 
did find a nonzero result for the pion scalar 
formfactor (ir\uu + dd\ir) at large q 2 in two- 
flavour HPChPT. Note that uu + dd is a sin- 
glet under SU(2)v but not under the full chi- 
ral SU(2)l x SU{2)r so the calculations given 
above do not restrict the chiral logarithms at 
large q 2 here. 



4 Comparison with experiment for 

XcO,Xc2 -» PP 

The fact that the chiral logarithms of order 
m 2 log m 2 vanish means that the leading term 
in the expansion of m 2 of the amplitudes van- 
ish. One could thus expect that the SU(3)v 
breaking in % c o,Xc2 — ► PP should be some- 
what smaller than the "usual" 20% as e.g. in 
Fk I 'F^ or mA/m p . Terms of order m 2 are how- 
ever not predicted, so a very clear statement 
that SU(3)v breaking effects are small is not 
possible. 

Let us however see how well the measured 
amplitudes live up to the "small SU(3)v break- 
ing effects." We take the input data from [T]. 
These are given in Table \T\ It should be taken 
into account that the 7T7t is the sum of 7r + 7r~ and 
7r°7T final state. Isospin predicts that the tt°tt° 



Johan Bijnens, Ilaria Jemos: Chiral Symmetry and Charmonium Decays to Two Pseudoscalars 5 







Xc2 


Mass 


3414.75 ± 0.31 MeV 


3556.20 ± 0.09 MeV 


Width 


10.4 ± 0.6 MeV 


1.97 ±0.11 MeV 


Final state 


10 3 BR 


10 10 G [MeV- 5/2 ] 


10 3 BR 


10 10 G 2 [MeV- 5/2 ] 


7T7T 


8.5 ±0.4 


3.15 ±0.07 


2.42 ±0.13 


3.04 ± 0.08 


K + K~ 


6.06 ±0.35 


3.45 ±0.10 


1.09 ± 0.08 


2.74 ±0.10 


K%K% 


3.15 ±0.18 


3.52 ±0.10 


0.58 ± 0.05 


2.83 ±0.12 


VV 


3.03 ±0.21 


2.48 ± 0.09 


0.59 ± 0.05 


2.06 ±0.09 


ft ft 


2.02 ±0.22 


2.43 ±0.13 


< 0.11 


< 1.2 



Table 1. The Experimental results for the decays Xrf>, Xc2 — > PP and the resulting factors corrected for 
the known m 2 effects. 



is half the 7r + 7r~ state and the data that go into 
the PDG average are compatible with this. We 
thus need to multiply it by 2/3 to get the tt + it~ 
result. Similarly, the KgKg final state is half of 
the K ~K° final state. Here we need to multiply 
by two to obtain the full final state. 

One question is how we deal with factors 
of m 2 P that we know are present. We include 
the phase space correction and note that for 
XcO the amplitude always contains p\ ■ pi — 
(m^, — 2m P )/2. The phase space contains the 

factor \pi \ = J m 2 , — Amp/2. For the decays of 
the scalar we thus define 

G = y/BR/\Pi\/(Pi-P2)- (8) 

The same question arises for the Xc2- Here the 
amplitude must always contain 'X^PinPiu- The 
field X2 must be replaced by its polarization 
tensor T// 11 and the amplitude squared and av- 
eraged over. The formula to perform such a sum 
is 

polarizations \ 

-^""/{""J (9) 

where K>* v = -g^ v ± q^q v /m 2 x and is the 
four-momentum of the Xc2 particle. The contri- 
bution thus always contains a factor 

pol 

(10) 

Alternatively we can choose the final state con- 
figuration with pi — (Ep,0,0,\p\\) and pi — 
(Ep, 0, 0, — |pi|) and choose an explicit set of 



five orthogonal polarization tensors satisfying 
n^T — f) — T — T and 

T x %Tx {h)llv = 5 ab which shows that the ampli- 
tude squared always contains a |pi| 4 . We thus 
define a normalized factor also for the Xc2 de- 
cays via 



G 2 = ^BR/\ Pl \/\ Pl \ 2 . (11) 

Looking at the columns Go and G2 for the 
decays to pions and kaons we see indeed that 
SU(3)y breaking is somewhat smaller than usual, 
about 10% for both XcO and Xc2 decays. For the 
decays to r\r\ we get a decent agreement in both 
cases but th^3 rfrf is not so good for the Xc2 
decay. 



5 Conclusions 

In this letter we have calculated the chiral log- 
aritms for the decays XcO, Xc2 to two light pseu- 
doscalars. We have found that they vanish and 
that this is a general result for all such decays 
of heavy chiral singlet states of spin and 2. 
We checked our result against the known result 
for the energy momentum tensor in ChPT and 
showed using the example of the scalar formfac- 
tors that the chiral logarithms do not vanish at 
high q 2 in all cases. 

Our result implies that the SU(3)v correc- 
tions are expected to be "small" in these de- 
cays and we have compared our results with the 
available experimental data and find that the 
corrections are reasonably small. 

3 The arguments given in this paper are only ap- 
plicable to the 7/ in the large N c limit where this 
state also becomes a pseudo Goldstone boson. 
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